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SYNOPSIS 


Uma Kapur 
Ph.D. 

Department of Physics 
Indian Institute of Technology, Kanpur 
August 1971 

CONFIGURATIONAL STATISTICS OF MACROMOLECULAR SYSTEMS 

Polymer chains consist of a large number of like 
units and because of the presence of internal rotations 
ziround single bonds^ are capable of taking almost an 
infinite number of conformations in space. 'This necessi- 
tates a statistical approach in studying their physical 
behaviour . The dimensions of polymer chains are most 
conveniently characterized in terms of the mean square 
end-to-end distsince ^ where n is the number of links. 

The knowledge of this distance is important in developing 
the equation of state for rubber elasticity, for deriving 
equations expressing intrinsic viscosity in terms of chain 
length, for interpreting light scattering data from dilute 
solutions and for promoting the growth of other important 
theoretical concepts. In the first chapter, we have 
discussed these aspects and briefly summarised the 
theoretical approach followed by other wcrkers in calcu- 
lating < R^* > . 
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In the second chapter, a reindom-walk model for 

hydrocarbon type chains with short range correlations is 

given. The three possible orientations of a hydrocarbon 

type chain have different probabilities depending on 

whether they are placed in a trans or in two possible 

gauche positions. The energy’ difference between a gauche 

2 

and a trans sequence of three links is taken as 6 . < > 

is calculated by describing a random walk on a diamond 

lattice. The progress of the chain is governed by two 

probability parameters. The mathematical methods of 

2 

Markoff chain are used to calculate < R^^ > for an oriented 
chain. For large values of n, following result is obtained. 


<R > 






^ 128 
4(2q+l-r)^ 


[3 ‘ [i+(<i-^)(i-q)3j 


It is shown that the above result reduces to those 
obtained earlier by Tobolsky and Eyring, The problem of 
excluded volume Is excluded. The Markoff chain method has 
been extended to semicrystalline polymers to calculate the 
chain dimensions. Generalisation of the ‘Tobolsky and 
Gupta' matrix has been done to calculate the expected square 
of the end-to-end distance, crystallinity, amorphicity and 
sequence lengths for semicrystalline polymers with orienta- 
tion in both the crystalline and amorphous regions. 

chapter III, the problem of phase transition 
between helix and random coil in polypeptide chains is 





discussed. Theory of Ising model has been used to study 
dimensional changes during melting in normal paraffines. 

It is based on the cell model for diffuse transition and 
involves writing the Ising matrix in terms of segment pair 
partition functions which are either in the amorphous or in 
the crystalline region. The experimental data of Flory and 
Vrij has been used to obtain the crystallinity vs. temperature 
curves. The nucleation parameter a is adjusted to give 
correct melting point . The dependence of transition curves 
on a values is also shown and, as expected, it gets sharper 
for smaller values. Density as a function of temperature is 
also calculated and experimental and theoretical results are 
compared. Because of the equivalence of the crystal struc- 
tures for polyethylene and N-paraffin hydroceirbon and the 
virtual identity of their unit cell dimensions, it is expec- 
ted that the limiting melting point fcr N-paraffin hydrocarbon 
series must be coincident with the melting point for linear 
polyethylene of infinite chain length. This is supported 
by the calculations. The nucleation parameter approaches 
the value 5 x 10~^, which yields correct results for poly- 
ethylene, Effect of pressure on the system is also studied. 
This Ising model formalism is extended to study the effect 
of pressure on. transition. The shift in melting point as a 
function of presstxre is evaluated and is in good agreeipent 
with the experinwntal results of Matsuoka. The theory of 



xi 


diffuse piiase transition has been applied to explain the 
thermal transition in three unit copolymeric system of 
L-Glutamic acid, L-lysine and L-alanine. Helicity vs. 
temperature curves are obtained at PH3 a^d PH8 values; they 
are also in good agreement with the experimental results of 
Morita, Simon and Blutt. 

The last chapter is devoted to concluding remarks 
and discusses various approximations made in the present 
Markoff chain approach and in treating the polymeric systems 
as a cooperative phenomena* 


I 
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INTRODUCTION 

Macr omolecular chains are formed by the successive addi- 
tion of monomeric units and have a threadlike structure. They 
consist of sequences of similar, though not necessarily iden- 
tical, chemical units connected together by covalent bonds 
into chains which frequently contain more than ten thousand 
atoms. They may exist as discrete molecules or they may be 
interconnected in various ways giving rise, in extreme cases, 
to a continuous network of chains in which the concept of indi- 
-vidual molecules loses its significance. Physical, chemical, 
thermodynamic and various other properties of these long chain 
molecules differ very much from those of the short chain and 
low molecular weight substances of similar composition. This 
is largely due to the size, shape and configuration that the 
polymer chain takes up. This, in turn, depends on the freedom 
of internal rotation around single bonds and the nature of side 
groups. The polymer may be either completely amorphous, semi- 
crystalline or highly crystalline according to the arrangement 
3rLd packing of the molecular chains. The mechanical behaviour 
of such systems depends on the distribution and extent of crys- 
"^^-llisation. In general, they can be classified as either 
natural or synthetic. In the former, we have, e.g. , cellulose 
3rid proteins whereas in the latter classification are included 
niaterials like polyethylene, polypeptides and polynucleotides. 

The subject matter of this thesis deals with a few 
problems regarding the configurational statistics, molecular 
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dimensions and the problem of phase transition in linear 

r\ 2 O 

hydrocarbons (finite system) and tricopolymers of glut" -lys - 
alan^^. The effect of pres-sure on the nature and transition 
in polyethylene is also studied. 

A polymer described by a chemical formula is a mixture 
of a multiplicity of different configurations. These, almost 
infinite number of conformations in space, come about because 
of the presence of internal rotations about single bonds. The 
problem of the geometry of the polymer chain, therefore, acqui- 
res a statistical character . The statistical treatment is 
based on the existence of large ensembles of members possessing 
independent degrees of freedom. Polymer chains are in fact 
such ensembles. However, the totality of polymer chains, -with 
the exception of a few, is never really an ensemble of identical 
objects. It is a complex mixture of entities which are hetero- 
geneous in composition and structure. Application of statis- 
tical theory to a given system generally requires a knowledge 
of its energy levels. But some polymer problems can be solved 
without this knowledge, e.g, , the state of a macromolecular 
system can be characterized by its dimensions and not by its 
total energy. Consequently, the states of the system can be 
regarded as states of different dimensions rather than states 
of different energy. 

In extremely dilute solutions of polymers, the indi- 
vidual chain molecules are virtually isolated from one 
another and the thermodynamic and gross physical properties 
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like viscosity, light scattering, etc. of the solution 

depend directly on the mean configuration and free energy 

of a single molecule. A quantity like the expected square 

2 

of the end-to-end distance <R > , where n is the number 

n ’ 

of links, of a polymer chain is important for determining 

the physical h^aviour of polymers. 

2 

The knowledge of > is important in determining 

light scatterdLng data function p( y ) ; which, as given by 
Debye ^ , is 


P(©) 



N 

I 


i,d*1 


Sin V R . . 

ii 

u R . . 


( 1 - 1 ) 


where © is the scattering angle, N is the number of atoms 
in the chain, R. j is the distance between ith and jth atoms, 

X vl 

y =: Sin and X’ is the wavelength of the light in 

the solution. Averaging is ceu^ried out over all chain confi- 
gurations.' By neglecting the higher order terms, following 
result is obtained 

P(©) = I exp (- y^ <R^j>^ /) (1-2) 

^ i,d 

The hydrodynamic properties of macromolecules in 
solution ^(viscosity , sedimentation, etc.) are determined 


by the dimensions of the macromolecules. For exan^lei, 
the viscosity of polymer solution is given by 



k 


/c2.3/2 

n = ^ ( 1-3) 

where M is the molecular weight and 0 is a constant. The 

2 

empirical value of 0 obtained by Flory by treating the 
results of experiments in which viscosity and light scattering 
were measured for different polymer solvent pairs, is equal 
to 


0 = 2.1 (+ 0 . 2 ) . 10 ^^ 

j” If (R^)^ is measured in cm. and n in cnP/g 

In this way, the above equation gives a convenient 

_2 

method of determining (R ) from the characteristic viscosity 

and the molecular weight . 

2 

The value of <R > is also useful in determining the 
coefficient of translational diffusion D. and the sedi- 
mentation constant S through the relations 



( 1 -^) 


M(1*»vP ' 

N.F. 
A t 


(1-5) 


(v is the partial specific volume of the polymer and p is 


the density of the solvent) , is the coefficient of 
translational friction and is related to the experimentally 


observed quantities. For linear chains in an ideal solvent 
(for which a gaussian distribution is valid) 
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= p < n ( 1-6) 

where p is a constant aind its value has been found to be 

2 

5.1. Hence knowing <R > and other quantities, the values 

of S and D, can be obtained . 
t 

2 

Earlier, the problem of finding <R > has been 
tackled by Eyring^, Lifson^, H^-eve^ and Suzuki^. They have 
assumed that the polymers are formed by the successive 
addition of C— C bonds in a random manner with the only 
restriction that the bond angles have specified values. 

An account of the calculation of chain length using the 

7 

method of matrix algebra has been given by Volkenstein , 

The properties of the polymers depend on the shapes 

or the conf igiirations that the polymer chains can adopt, the 

flexibility of the chains, and the presence and extent of 

the cross-linking or branching. The idea of the "excluded 

8 

volume effect" was suggested for the first time by Flory 

in 19^9* The excluded volume arises due to the long range 

9 

correlations of the chain elements. Montroll , Hermans 

10 11 12 13 

et al , Stockmayer and Fixman , James , Saito and 

1 4 

others have given due consideration to this * effect^ in 
developing their theories. Lifson ^ , Nagai and Tobolsky 
have considered the effect of branching on polymer dimen- 
sions. It is intuitively obvious that if the molecule 
consists of branches, the molecule is going to be more 
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compact than a straight chain of the same molecular weight . 
One can define a parameter g as the ratio between the average 
square length of a branched molecule and the average square 
length of a linear molecule . 


g 



> 


> 


( 1 - 7 ) 


o 

In general, g is less than one if there is branching 

1 8 

in the molecule. Lifson also considered the effect of 
charged side groups and determined the molecular dimensions 
for polyelectrolyte molecules. 

19 20 

Since the introduction by Kuhn ’ of the random 
walk model to approximate mathematically the distributions 
of the configurations that flexible linear molecules might 
assume, the theoretical calculations of the mean dimensions 
of polymeric molecules has been the subject of numerous 
publications^ ' ^ ^ ^ in recent years. In a random walk 
model, the actual molecule is replaced by an artificial one 
in which instead of having the molecule composed of actual 
bonds with fixed valence angles and a certain amount of 
restricted rotation, the molecule is constructed of statis- 
tical segments. It is assumed that a polymer chain is formed 
by the addition of successive bond vectors in such a way as 
if they are the successive steps of a random walk. The 
formalism of the mathematical method of Markoff chains 
has been used. Its use has been earlier demonstrated by 
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several other workers, e.g. , Montroll , Tobolsky and 
22 

Gupta . A polymer chain can be considered as a Markoff 
chain because the relative orientation of a unit is 
dependent upon the position of the preceding one. 

In the next chapter, two models based on random walks 
on cubic and diamond lattices are considered. In the former, 
a polymer chain with orientation in both crystalline and 
amorphous phases is studied whereas the diamond lattice model 
is used to approximate an oriented hydrocarbon type chain with 
short range correlations. A polymeric system is characterized 
by parameters such as expected squaire of the end-to-end dis- 
tance, the extent of crystallisation and crystallite sizes. 

These parameters have been evaluated. Of the two models, 
diamond lattice model is a better approximation in the sense 
that in linear polymeric chains such as hydrocarbons the 
valence C!— C bond angles are the same as in a diamond lattice. 

The third chapter is devoted to a study of order- 
disorder transition in macr omol ecul ar systems. Transitions in 
finite and infinite systems such as n-paraffins and polyethy- 
lene respectively are derived. In the latter case, the effect 
of pressure on transition is discussed. The Ising model approach 
is used and is modified to explain helix-coil transitions 
in a tricopolymeric system also. Order -disorder phenomenon 
and accompanying phase transitions have attracted physicists 
for a long time. They are of special interest in the case 
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of biopolymeric chains because the translation and trans- 
cription of the genetic information is closely related to 

the uncoiling of a DNA molecule. In tricopolymer ic systems, 

2 5 

it was observed experimentally- by Doty et al that some 
physical properties like the intrinsic viscosity and optical 
rotation, etc. underwent considerable change on changing the 
conditions like the temperature and solvent composition. 

This change in the properties due to change in such physical 
conditions was correlated to the change from helical structure 
to a random coil struct^u^e . This phenomenon is known as the 
"helix-coil transition" . As the transition is quite sharp, 
it is referred to as a phase transition of the first order 
because of the accompanied changes in the entropy of the 
system. 

Several attempts have been made to theoretically 

interpret this helix-coil transition and they are all based 

on a modification of the Ising model proposed for ferro- 

26 

magnetism. It has been studied by Zimm and Bragg , 

27 28 29 

Gibbs and Di Marzio ‘ , Rice and ¥ada , Lifson and Roig 

30 3 1 

Hill and others , Just as a lineair Ising chain is made 
up of two spin states defined as ‘up* euid 'down* , so is a 
polymeric system made up of 'helix* and *coil* states or 
'ordered* and 'disordered* states. The aim is to evaluate 
the total partition function for the entire chain in terms 
of contributions made by individual segments. Other thermo- 
dynamic properties can be obtained from a knowledge of the 
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grand partition function. It has been shown by Landau and 

32 

Lifshitz that a one -dimensional system cannot exhibit 
phase transition; however, experimentally it is seen that 
polymeric systems do exhibit what may be called a phase 
transition. The contradiction arises because of the assun^)- 
tions made by the above authors in their theoretical approach. 
The assumptions axe 

(a) the interaction has a finite range, and 

(b) the interaction is non-singular, i.e., of finite strength. 

Once these constraints are lifted, the one -dimensional 

system shows all characteristics of a phase transition. The 

cooperative nature of hydrogen bonding in biopolymeric systems 

takes care of the first restriction and the introduction of 

a nucleation parameter which is a measure of difficulty of 

transition from the random coil to the helical phase in a 

way introduces singularity in the nearest neighbour interaction. 

The sharpness of the treinsition is governed by the magnitude 

of this parameter. It becomes more and more step like for 

33 

decreasing values , Baker^^ has also shown that a one- 
dimensional system will approach even a second order transi- 
tion in the limit tiiat the range of its interaction becomes 
infinite. The nature of the transition also depends on the 
presence of the defects aind size and distribution of ordered 
regions. The width of transition for finite system is, in 
general, more than that for infinite ones. 
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N-paraffina with finite degree of polymerisation and 
linear polyethylene with almost infinite chain length belong 
to a different class of polymers. Eventhough it is difficult 
in these cases to define cooper ativity factors as precisely 
as in hydrogen bonded biopolymers, one could still use the 
formalism of I sing model . The heat of transition in such 
cases would reflect both inibermolecular and intramolecular 
interactions and this in a way also takes care of the three- 
dimensional structure eventhough the model used is a one- 
dimensional one. Because of the equivalence of the crystal 
structures for polyethylene and N-paraffin hydrocarbon and 
the virtual identity of their unit cell dimensions, it is 
expected that the limiting melting point for N-paraffin 
hydrocarbon series must be coincident with the melting 
point for linear polyethylene of infinite chain length. 

This is supported by the calculations. The nucleation 
parameter approaches a value which is approximately the 
same as for polyethylene. This problem has also been 
studied by Flory and Vrij^^, Broadhur st^'^ and others^^. 

Flory' s treatment is based on the supposition that both the 
enthalpy and the entropy of fusion may be adequately repre- 
sented by linear expressions, each comprising a term 
proportional to n and an additive constant for contributions 
associated with the end groups. The melting point T^^^ is 
expected on this basis to be related to the chain length 



n as 


1 1 


= AH(n+a)/ AS(n+-b) = (n+a)/(n+b) ( t-S) 

where AH and As represent, respectively, the enthalpy 

and entropy of melting per CH2 group in the limit n -»■ <» , 

and aAH and b As are the end-group contributions. It is 

to be noted, however, that and other parameters appearing 

in equation (I-8) as well, will depend, in general, on 

temperature and should not, on the basis set forth above 

for this equation, be treated as constants. Revisions 

appropriately taking account of the dependence of on 

35 

ten^jerature were introduced by Broadhurst , His analysis 
retains, however, the basic assucption of linearity of the 
enthalpy and of the entropy with n at constant temperature. 

Lastly, a part of Chapter III is devoted to a study of 
pressure. ^e transition curve . The mechanical behaviotar 
and morphological structures of a polymeric material are 
markedly changed when subjected to stress, e.g, , vulcanised 
rubber which does not crystallize at ordinary temperature 
becomes highly crystalline when subjected to tension. 

Further the elongation of a network polymer in the amor- 
phous state increases its free energy and if the polymer 
is capable of crystallising, its melting point will also be 
increased. The effect of strain on melting point has been 
theoretically studied by Flory^"^. His model consists of 
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bundle like crystallites with the direction of the chain 
through the crystallite being perfectly aligned along the 
drawing direction. Furthermore, the exit from one crysta- 
llite, and the entrance into the next, was assumed to occur 
in the raainner which minimized the path of the intervening 
amorphous chain. These conditions imply that the sanple 
is rather well drawn prior to crystallization. The equi- 
librium degree of crystallinity, o) , for any ten^erature T 
was found by letting ( SF/ )_ vanish, and the melting 
temperature T^ was obtained from this relation by setting 
05 =. O . This model is somewhat oversin^Jlif ied because it 
fails for small strains as the assumption of peffect 
alignment is no longer valid. In the present approach, the 
formalism of Ising model in a magnetic field has been used. 
The relative contribution of a segment in a crystalline 
state with respect to one in the amorphous state is changed 
due to applied pressure and this is built into the Ising 
matrix. The theory of cooperative phenomenon is used. 

The last chapter is the concluding part of the 
thesis. Here, a brief comparison of various approaches, 
limitations of the underlying assun^jtions and their effect 
on the results obtained are discussed. Suggestions for 
further work and in5)rovement are made. 
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CHAPTER - II 

MOLECULAR DIMENSIONS IN SEMICRYSTALLIKE AND 
ORIENTED CHAINS 

2.1 INTRODUCTION 

As discussed briefly in the first chapter, one of the 

parameters characterising a macr omolecular chain is the 

2 

expected square of the end-to-end distance <R^ > , where n 

is the number of links. This distance is important in 
developing the equation of state for rubber elasticity, for 
deriving equations expressing intrinsic viscosity in terms of 
the chain length, for interpreting light scattering data from 
dilute solutions and for promoting the growth of other impor- 
tant theoretical concepts. Various methods have been used 
for its calculation depending on the geometry of the molecule. 

In this chapter, Markoff chain method has been developed 

to derive explicit formulas for the n step transition proba- 
Ic 

bilities p which are used in the calculations of the 
mn 

expected square of the end-to-end distance for a polymer 
chain . 


2.2 MARKOFF CHAIN METHOD 

In the theory of Markoff chains, the simplest generali- 
zation is that the outcome of any event depends on the outcome 
of the directly preceding event (and only on it). In the 
present problem of the calculation of molecular dimensions 
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of the polymer chains, the above-mentioned generalization 

has been imposed on the random walk model so that Markoff 

chain method can be applied to it. As comprehensive reviews 

1 2 

of this method are available in literature ’ , only a brief 

account of this method is presented here for the sake of 
completeness . 

Let a Markoff chain be designated in terms of a set of 
dependent variables Z^, Z^ , Z^ and let these variables 

be related so that the probability distribution of Z^ 

(j = 0,1,2,..,) is uniquely defined by that of ^ and the 
law of dependence is independent of j. 

Let the possible values of Z^ be a^ , a.^, , . . Then 

the probability distribution of Z . is defined by 

cl 

= pr(Zj = aj m= 1 , 2 ,...N 

and j = 0,1,2,... ( 2 - 1 ) 

(Here the common notations of the theory of proba- 
bility have been used and the symbol pr(Zj = a^) denotes 
the probability that Z^ equals 

As Z . can take one of the values from the set 
3 

a^ , ag , . . . we know 

I ~ ^ ^ (2-2) 

mssi 

A single Markoff chain can be completely characterized 
in terms of the initial probability distribution of Z^, i,e., 
p^(o) and the transition probabilities p ^ 
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p (n m) 




if Z . 
3 


aj (2-3) 


p > o xf 1 < n, m < N 


Since Z . = a implies ^vith certainty that Z . . has 
.1 n ^ , 14-1 


one of the values a^ ,a^ , * . .a.., we have 

1 ' 2 N 


N 

I = 1 for n = 1,2,...N (2-4) 


The probability of the simultaneous occurrence of two 
linked events is equal to the product of the unconditional 
probability of the occurrence of one of the events and 
conditional probability of the other assuming that the first 
event has occurred. Hence 


pr(Z = a and Z. = a ) = p P„(j) (2-5) 

Since the probability of occurrence of a given event 
is the sum of the probabilities of mutually independent ways 
in which it can happen, 

N 

Pni(j + ’') = I Pmn ^ ^ 

n=1 

... ( 2 - 6 ) 

Similarly the distribution of is obtained by 

iteration of (2-6), as follows 

N 

P„(J+2) = I PmiPn<j+’) 

n=1 


(2-7) 
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In general 


N 

= I 

n^=1 

11 

« • » « 

N 

Ip p 

n^=1 ”l”2 

m s 0 0 







. . .(2-8) 

Hence can be 

calculated from 

initial distribu- 

tion Pjjj(o) [when j=o 

] by 

application of equations (2-7) 

and (2-8) . 






This multiple 

summation process can 

be avoided by 

writing* these relations in 

a matrix Term. 

The transition 

probabilities p can 

•*^mn 

thus 

be arranged in i 

a matrix of transi- 

tion probabilities, 






P11 

P 12 

P*!^ • . . # 

^1N 



P 21 

P 22 

P 23 • • • ‘ 



p = 

P 31 

• 

• 

P 32 

^33 • • • • 

i 

! 

(2-9) 


# 

Pni 

%2 

Pn3 • • * • 

%N 



where the first subscript stands for row and the second for 
column. 


P is a square matrix which has non-negative elements 
and the sum of all the elements in a row is unity. Such a 
matrix (finite or infinite) is called a stochastic matrix. 
Any stochastic matrix can be used as a matrix of transition 
probabilities ani with initial distribution { Pm(o) } it 
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completely defines a Markoff chain. 

, 2^ be the set of row vectors. 


Let X-i > 


• * « • 


, . . . . , be the set of column vectors, , . . . . , 

be the corresponding- eigenvalues of the matrix P. These 
quantities satisfy the equations 


^ ^ ^ ik = ^k ^k (2-10) 

k 

is the kth step transition probability from the 
state m to n. These elements can be evaluated from the 
following relation 



N 

I 


r = 1 




( 2 - 11 ) 


C“ is the inverse of the inner product of row and 
column vectors. 


2,3 EXPECTED SgiJARE OF THE END-TO-END DISTANCE 

Consider the freely rotating segments of a chain to 
be each resolved into three components of length ^ along 
X, Y, Z axes. These are called link vectors. Let a chain 
of n links be represented by vectors ( ), 

i = 0,1,2,,., n-1 , The expected square of the end-to-end 
distance is denoted by 

<R^> = E [ ( x(°) + X^^^ + + )^ 

+ ( + Y^^^ + + )^ 

+ ( z(°) + Z^^^ + .... + )^] 


( 2 - 12 ) 
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3n+2 

n-.2 

I 

n-1 / . 

1 E [X^^^ 



i=o 

j=i + 1 


3n+2 

n-2 

1 

n-1-i / V 

I E [X^°) 



1 = 0 

k=1 



Substituting 

= E + Z^°) ^(2-l4) 


We get , 


<R^> 

n 


n-2 n-1-i 

= 3n + 2 I I L 

i=o k=1 
n-1 n-1-i 

= -3n + 2 J I Lj 

i=o k=o 


(2-15) 


There can be large number of configurations which a 
molecular system can assume. Therefore, one has to apply a 
statistical approach. It is convenient to replace the actual 
molecule with a synthetic one in which the actual chemical 
bonds are replaced by artificial ones. Instead of having 
the molecule composed of actual bonds with fixed valence 
angles and a certain amount of restricted rotation, the 
molecule is constructed of statistical segments. These aire 
chosen such that the average square length of synthetic 
molecule is the same as the average square length of the 
real molecule. For this pxirpose, models are chosen which 
are nearest to the stmcture of actual molecular system. 
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To calculate the moleculajr dimensions of polymer chains, 
two types of models are considered. In the first, the polymer 
chain traces cut a path on a cubic lattice (Fi^re l) and in 
the second, the polymer chain traces out a path on a diamond 
lattice (Figure 2). In the former, the general problem of 
molecular dimensions in a semi-crystalline polymer chain with 
orientation in both ordered and disordered phases is consi- 
dered, whereas in the latter, the conf crmational statistics 
of an oriented hydrocarbon type chain with short range inter- 
actions is studied. 

2 *4 CUBIC LATTICE MODEL ¥ITH ORIENTATION IN AMORPHOUS 

AND CRYSTALLINE PHASES 

Real polymer chains consist of amorphous and crysta- 
lline regions * The sequences of segments in these regions 
can be represented as 

r r r hhr hhrhr hhr * 

where r corresponds to a segment in the amorphous region 
and h that in the crystalline region* L^t us assume that 
the semicrystalline polymer chain under consideration traces 
cut a path on a cubic crystal lattice such that the crystalline 
and amorphous regions of the chain have preferred directions 
only along +x, +y and axes. Each segment consists of 
three links pointing successively in the three preferred 
directions. 

A +x crystalline region is one in which the chain 
progresses in the 4 -x direction but in which "the y and z 






■^f0,2.HYI^OCARBON TYPE CHAIN ON DIAMOND LATTICE 



24 


vectors are symmetrically spiralled so that the net contri- 
bution of the y and z chain components over this entire 
region is zero. Similar discussions are applicable for crys- 
tals in the other five possible directions. If the bond 
length or the length of a segment is taken to be l^j then the 
length of each link in the amorphous or crystalline regions 
will be 1 q / 1/3. But for simplicity this length is taken to be 
unity in calculations and the factor l^/v^ is again introduced 
at the end of calculations . 

The succession of ±x links (i.e., every third link) 

of the polymer chain forms a Markoff chain. Consider the 

elements of a chain to have five possible states r^, r , 

h , h and h . The state r denotes an x link in the 
+ - o + 

amorphous region pointing in the +x direction, r denotes 
an X link in the amorphous region pointing in the -x direction, 
h^ denotes an x link in a +x crystalline region, h denotes 
an X link in a -x crystalline region, h^ denotes an x link 
in a ±y or ±z crystalline region (the average length of such 
a link is zero) . 

Tobolsky and Gupta have treated the sequence of x 
links by Markoff chain method taking into account the 
orientation only in the crystalline region. They have used 
the following transition probability matrix for describing 
the succession of states. 



25 



r 

r 

h 

h 

h. 


+ 


-t- 


o 

r 

+ 

Ki-p) 

i(i-p) 

^x ^ 

0 

(1- Up 

r 

i_p) 

i(i-p) 

0 

^ P 

(1- B^)p 

h 

1 - a 

0 

a 

0 

0 

+ 






h 

0 

1 - a 

0 

a 

0 

h 

o j 

^(1 - a ) 

1(1- «) 

0 

0 

a 


(2-16) 


■where the parameter p represents the probability that a 
bond in an amorphous region will be followed by a bond in 
crystalline region, the quantity a is the probability that 
a bond in a crystalline region is followed by another bond 
in the same crystalline region and the factor is the 

fraction of crystallites oriented in the +x direction. 
Similarly and are the fraction of crystallites 

oriented in the +y and +z direction respectively. Therefore, 


8^ + By * 6^ = 1 (2-17) 

For an unoriented crystalline polymer, 

8^ = By = % = 1/3 (2-18) 

and for complete orientation in the x direction, 

= 1, = 0 (2-19) 

Tobolsky and Gupta^ have further calculated the expected 
square of the end-^to-end distance projected along the x, y 
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and z directions, average sequence length in the crystalline 
and amorphous regions, fractional crystallinity and fractional 
amorphicity using the Markoff chain method. All these results 
were derived for a chain of n links in each of the three 
directions, i.e,, for a chain of n bonds » These results are 
given below* The expected square of the end-to-end distance 


in the x direction, <R > for n links 

X 


<R-> = n 


1+P ) 

TTTTp" 1-p g„ (1-p - “ )C - “+p] 


1-a + 


-rl 


o 

3 


. . . ( 2 - 20 ) 


Similar expressions with 3 and 3 can he written for 
^ ^y z 

2 2 
<R > and <R > , 
y z 

The sequence length in the amorphous region 
(uninterrupted by r^ or r states) 


= ■•/(■<- a ) 


( 2 - 21 ) 


The sequence length in the amorphous region (uninterrupted 
by h or h states) 

= 1/p (2-22) 

Fractional crystallinity = p/(‘'-a+p) (2-23) 

Fractional amorphicity = ( 1 - a )/( 1 - a +p) (2-24) 

Here the orientation in the amorphous region is not 
considered. However, there are polymers which show orienta- 
tion in amorphous region also as has been experimentally 
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observed in X-ray studies^ in samples sucb as polyethylene^ 

3 

In the present work Tcbolsky and Gupta^s work has 

been further extended by including the orientation in the 

\ 

amorphous region also. Now the transition probability- 
matrix is given as 

r r h h h 

+ _ + - o 

r^ y 1-p-V a p 0 (1- B )p 

r_ (1-p-p ) V 0 3..P (1- 6 )p 

p-h^ (l-a) 0 a O 0 ( 2 - 25 ) 

h_^ 0 (1-a) 0 “ 0 

ho |_(1- a)/2 (1- a)/2 0 0 a _ 

Here an extra parameter ja is introduced, which denotes 
the orientation in the amorphous region. To calculate the 
molecular dimensima Markoff chain method has been used and 
only those components are considered which are in the -fex 
direction. If the length of the successive vectors be 

x^^^ then the expected square of the 

end-to-end distance is given by 


^ 2 
< R > 

X 

r (o) 
= [ x'- ' 


(n-1 ) 

• + x'' ^3 

( 2 - 26 ) 


n-1 , . ^ 

= I x<") 


n-2 n-1 

. 2 I 1 

^(i) Ji) 


isso 

n-1 

ip:o j=i+1 

. . .( 2 - 27 ) 


= nE(tl) + 

2 1 
k=1 

(n-k) 

(2-28) 



where E(± 1 ) is the expectation that an x vector has the 
value 11 rather than zero and Lj^ = E( x^°^ x ) . 

In the Markoff formalism the value of L, can be 

xC 

written as 


= Pn - Pi2 + p'fa - (2-29) 

Id k. 

where ^ , P■^ 2 ’ etc., arc the elements in the first row of 
the kth power of P and it is assumed that the chain started 
in the r state. 

* 4 " 

The chain can start in any of the five possible 
states, i.e., r^, r^, h^ , h^ or h^ . The growth of the chain 
will depend on its initial starting state. The end-to-end 
distance is calculated for the chain starting in each one of 
the frve states separately and are weighted according to the 
respective steady state probabilities. They are finally 
combined to get > . Using the notation 1,2, 3, 4 , 5 for 

s^nd h^ respectively, we have 



2,5 STATIONiffiY STATE PROBABTLIIY 


(2-30) 


¥hen N is sufficiently large the system eventueilly 
attains a steady state and there is a finite probability 
called stationary state probability, of its being found 
in any one of the available states irrespective of the 
initial state. In statistical mechanics this corresponds 
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to the system which has attained equilibrium; i.e., there 
is no xurther time development of the system. These sta- 
tionary state probabilities would then be the normalised 
components of the eigenvector corresponding to eigenvalue, 
X = 1. 


The stationary state probabilities p are given by 


the relations = m and the stochastic property ^ ^ij ~ ^ ' 

1 - a 


2 ( 1 - a +p) 
BP 

2 ( 1 - a +p) 
P(1- B, ) 
(l- a +P) 


f or 

r 

+ 

and r states 

(2. 

- 31 a) 

for 

h 

+ 

and h states 

(2- 

- 31 b) 

for 

h 

o 

state 

(2. 

- 31 c) 


From these 


■ E(±1) = 


(l-a + 3p) 

1 - a +P 


(2-32) 


Hence from equation (2-28) 

1 -a 


<R^ > 


2(1 - a +p) 


1 - a + Up n-1 

n :: + 2 I (n^k) L, 

k=o 


1- a+P 


(2-33) 


Similar expressions can be written for the other 
states. The problem thus reduces to the evaluation of Lj^. 
which involves the kth-step transition probabilities which, 

t 

in turn, can be obtained with the help of equation (2-11). 

The matrix of transition probabilities has eigen- 
values 1, a , a-p, X' and X", where V and X" are 
given by 
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X' + 

X" - 

a + 2y - 1 + p 





X' 

x" 

-p ( P - a - ag ) 

+ a (2y 

-■>) 



Ic 

The various terms oT p are 

^mn 




k 

^11 = 

k 

1 

A X^+'' . 

. a A X"'^ 


(2 -34a) 

^22 

2(l-a+p) 

2 AX 


k 

k 

1 

A X^+'' - 

A 

- a A X 


(2 -34b) 

^12 " 

^21 

2 ( 1 - a +p) 

2 AX 


k 

k 

P 24 

Spl- 

8P_ 

- A X^ 


(2-34c) 

^13 ■" 

2 ( 1 - a +p) 

2A X 



k 

■r» . 

k 

gpl' 

.__B_P 

- A X^ 


(2-34d) 

Pl4 - 

P 23 

2(l-a+p) 

2 A X 

“ ii A 


k 

^15 ~ 

k 

P 25 

(1- 6 )pl’ 

(1- a+p) 




(2-34e) 

k 

P 31 = 

k 

P 42 

(1- a U’ 

- 2 ( 1 -a+p) + 

( 1 - a) 
2A X 

ax'" 


(2-34^) 

fU 

(1 

k 

P 41 

1 I- “ ) v 

2 ( 1 - a +p ) 

(1- a) 
24 X 

ax'" 


(2-3^g) 


k 

P 44 

g 1" 

11- BJ 

, a A X^- 

ak _ . _ 

^4 

, . ,k-1 

X^ A X 

^33 " 

- 2(1- a+p) 

2 

a + ^ 

A X 

. . .(2-34h) 


4^ 

ii 

k 

31" 

(1- 3) 

a A Xk 

nk 


X AXk-1 

5 

P 43 

2(l-a+p) 

2 

“ - 2 

A X 

. . .(2-34i) 


where 

1 

= (1- a) + p(a 

-P) 





1’ 

= 1 -(a -P)^ 






1" 

= p + (1- a ) (a -p)^ 




and A 

is the 

operator for X 

’ - X" . 





and 
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Substituting equations (2-3^a - d) in equation (2-29) > 


L,j^ = 5 (t> + (BP-a) D ) (2-35) 


Making use of the following equations, 


n-1 

I (n-k) X 

k=1 


1 - X 


nX X^^+^ 

(1- X ) + ^^2 ^^2 

dx 

( 1 - a) (2-2 y -p- g p) 


(2-36) 

(2-37) 


D _ D>.(e^^.2P-Up'ltp(B-c^- »e^-a(2u-1> /, 

1-X - (1-0 ) (2-2y - Bp) - P ^ 


and retaining only the leading terms, one gets 

n-1 n [(1- a )(2 V +p.l) +Bp(2-(x)] 

I (n.k)L = 

k=1 ^ (1-a )(2-2y - p - 6p ) 


(2-39) 


In the manner outlined above, it can be easily seen 


that 






<r 2 > 

*2 

= <R^ > 




( 2 - 40 ) 

<R^ > 

"^3 

= <R^ > = 

Bpn 


- ( 1 - a+ Bp ) 


2 ( 1 - a p ) 


_ ( 1 - q + p) 






2 [ 1 + Bp ( 1 - a)- 

o (2 V - 1 +p )] 





( 1 -a )( 2 - 2 y -p 

- %) ) 






...( 2 - 41 ) 

<R^ > 

^5 


1 - a + Bp 

n 


( 2 - 42 ) 

( 1 - a +p) 

1 - a + p 
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Thus, finally, from equation (2-30) 


<R > = n 

X 


+ S'!) (2-2 y -p + gp ) 


2 ^ (2- a) 


.1-a+ p (2-2]i -p - gp ) ( 1 -a) ( 1 -a+p) (2-2y-p-3p) 


2(p +2y -1) [(1- a)' 


a3p] 


<R^ > 
z 


( 1 - a H-p) (2-2y - p - 3p)(l-a)J 

...(2^43) 

2 

Similar expressions can be written for < R > and 

y 


2.6 FRACTIONia. CRYSTALLINITY AND AMORPHIC ITY 

The fraction of x components which exist in the h^ , 
h and h^ states can be regarded as the fractional degree of 
crystallinity and is given by 


(P -13 + P-14 + P 15 ) - +p) (2-44) 

Similarly the fractional degree of amorphicity is given by 

l^™k->oo(Pii + P- 12 ) = (i _~a +p) (2-^5) 


2 .7 SEgUENCE lengths IN CRYSTALLINE AND AMORPHOUS REGIONS 

The probability that a vector continues in the 

n 1 

cxystalline + or - state over (n-1) steps is a ” and that 
it gets out in the nth step is (1- a). Hence sequence length 
in the crystalline state is 

lim 2 n(l-a ) 0 ”“' = / ^ \ (2-46) 

n-»-oo ' ' ^ 1 _a) ' 
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SimilEirly, sequence length in the amorphous region (uninte- 
rrupted by crystalline states) 

= 1/p (2-47) 

By substituting p = , the results ot Tobolsky and 

3 

Gupta can be obtained. 

2 .8 DIAMOND LATTICE MODEL FOR A HYDROCARBON TYPE CHAIN 
WITH SHORT-RANGE CORRELATIONS 

Now consider that the polymer chain traces out a path 
on a diamond lattice. This lattice is a better approximation 
than a cubic lattice in the sense that in some polymeric 
chains such as linear hydrocarbon t 3 rpe, the valence C— C bond 
angles are the same as in a diamond lattice. Various polymers 
have such bonds, for example a polymer like polyethylene has 

5 

isomers shown in figure 3. In a diamond lattice each 
lattice point is joined to its four nearest neighbours 
through tetrahedral bonds (Figure 4a) and, therefore, each 
carbon atom has four such vectors ending on it. The aim 
here is to find a set of axes which can describe these bonds 
in a convenient way. For this an orthogonal set of axes is 
chosen and the points A, B, C, D are defined by (0,0,0), 
(2,2,0), (0,2,2), (2,0,2). Joining these points gives a 
tetrahedron thereby allowing us to label A, B, C, D as atoms 
(Figure 4a) . 



X 





36 


Th.e coordinates of the center (where the perpendi- 
culars from any vertex to the opposite face meet) will he 
given by 

/ 0 + 2 + 0+2 , 0 + 2 + 2+0 , 0 + 0 + 2+2 i 

j. _ _ j or I , I , I j 

Hence the vector AE is given by (1,1,1). Similarly 
other vectors are EB (1,1,1), EC (l,1,l) and ED (1,1,1 )• 

In the case of an actual lattice, other vectors are also 
present but they are described with respect to the axes shown 
in Figure 4b . They will be given by (1,1,1), (1,1,1), (1,1,1) 
and (1,1,1). Thus, a set of eight vectors is obtained, which 
describe all the possible bonds. In solving the problem by 
the method of ’random-walk' , these vectors also describe the 
direction of translation in going from one lattice point 
(atom) to the neighbouring one. This set can be symbolically 
written as (1,1,1). From here, a property of the diamond 
lattice emerges, that each successive vector is obtained 
from the preceding one by changing the sign of one component 
amd leaving the two others unchanged (a property which makes 
the Markoff chain method applicable to the system). The 
three possible orientations of a link have probabilities a 
or b according to whether it is placed in the trans or in 
one of the two gauche positions with respect to the two 


preceding links. 
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The probabilities a and b are determined by the 
relations 


a + 2b = 1 (2-48) 

f = (|t) (2-49) 

where e = E -E, is the energy difference between a gauche 
g rr 

and a trans sequence of three links. 

Earlier Tobol sky^ calculated the expected square of 
the length of a hydrocarbon type chain by considering the 
random walk on a diamond lattice. He considered the proba- 
bility for two successive changes in sign to occur at the 
same coordinate [i.e., is identical 

to (X^, Y^, Z^) ] to be r and the probability for the second 
change of sign to occur at one of the other two places with 
probability (l-r)/2. These three situations correspond to 
the trans and the two gauche configurations. 

Now substitute 




= (1,1,1) 

(2-50) 

and 

u. = . y(^) 

+ 

(2-51) 


Then from equation (2-l4) 

= E(U^) (2-52) 

The variable U. can assume the four values +3i +1. 
Tobolsky considered the sequence of pairs (tJ^#U^) , (U^jUg). 
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(U^jU^) (U^ 1 each padLr 

corresponds to a state with possible values (3 1), (l 3), 
(1 _i ) (_1 1) (_i _ 3 ) and (-3 -1). The transition proba- 
bility matrix can be given as 



The expected square of the end.-to--end. distance- 
for large values of n is obtained as follows, 

<R^> = (2+6r) nl^ /3(l-r) (2-54) 

In the present work a further modification of the 
Tobol sky ' s model is considered. A sequence of three 
consecutive links is taken. If the last one is parallel to 
the first, i.e., when the sign of the same coordinate has 
changed, it is assumed to occur with a probability r. This 
is the probability of the orientation being destroyed. Let 
the other two probabilities-, i.e., of preserving the orienta- 
tion (in the direction indicated by the arrow in Figure 2) 
and for reversing the orientation be designated as p and q 
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respectively. In other words, if two consecutive changes 
of sign involve the coordinates X and Y, then the next change 
occurs at the coordinate Y again with probability r, at the 
coordinate X with probability p and with probability q at 
the remaining coordinate. 

Now the transition probability matrix is given below 


3 1 

1 3 

1 -1 
-1 1 
-1 -3 

-3 -1 


3 1 
0 
1 

0 

0 

0 

0 


1 3 
r 
0 
0 

q 

0 

0 


1 -1 
1 -r 
0 
0 

p4r 

0 

O 


-1 1 
0 
0 

p+r 

O 

0 

1-r 


■1 -3 
0 
0 
q 
o 
o 


-3 -1 
0 
o 
0 
0 
. 1 
0 


(2-55) 


Obviously p + q + r = 1 . 

From equation (2-15) expected square of the end-to-end 
distance is given by 


<R^ > 
n 


where 


n-1 

= -3n + 2 I (n-k) 

k=:0 


(2-56) 


The chain can start in any one of the six states. 
The weights assigned to the various states etre 



ko 


03 


1 


3 


03 


2 


03 


3 


03^ = 


03_ 

5 

“6 


= 1 
= -1 
= -3 


With these , in general 

^(^(s) ^(t) ^(t) ^t)) 


V ( s) I 

. I . Pij “j Pf 


(t-s) 


d.k=1 


3 


( 0 , 


...(2-57) 

where is the step transition probability from 

-I- J 

the state i to j. The matrix P has three eigenvectors of the 
symmetric pattern (x,Y,Z ,Z,Y,X) > which will not contribute 
to Lj^ and three eigenvectors of the antisymmetric form 
(X,Y,Z , -Z , -Y, -X) which would contribute to Lj^. The eigen- 
values corresponding to the ant isynune trie eigenvectors are 
- 1 , and , wher e aiid 


X 2 are roots of 


X-| “3" X^ “ ' Q. 
Xi X 2 = 


The p(^^ elements which contribute to L, are given below 

X J| K 


p(k) _ 
p-l 

(k) 

“Pl6 

(k) 

= -P6l = 

(k) 

P66 


(2-58a) 

k 

^12 

k 

-Pi 5 

k 

= -P62 = 

k 

P 65 

= (- 1 )^+^ C,+ C 2 + 

Xg'^^C^ ( 2 - 58 b) 

k 

k 

k 

k 

/ ,^k +1 1 -r .k,. 

- (“U — ^ 1 + ^1 (1 

“ ^1 ) C 2 

^13 “ 

-Pl4 

= -P63 = 

P64 





+ XgCl- ^ 2 )^ 

...(2-58c) 



4 l 


k 

k 

k 

k 

k 

P21 ~ 

-P26 = 

-P51 = 

P51 

- P56 ““ 

k 

k 

k 

k 

k 

P22 ■" 

P1I = 

“P25 “ 

■‘P52 

- P55 

k 

k 

k 

k 

II 

+ 

P23 = 

-P24 = 

-P53 - 

P54 

k 

k 

k 

k 

- (-1)'^- 

+ 

P3I - 

-P36 = 

*-P 4 l = 

P46 

k 

k 

k 

k 

k 

P32 “ 

-Pi 3 - 

-P33 = 

“P42 

= P45 

k 

k 

k 

k 

= (-1)’^ 

P33 - 

-P34 = 

1 

ii 

P44 


,k+i 


.k-1 


k-1 

2 ^3 


. . .(2-58ci) 
(2-58e) 


Ci+ (1- 


(2-58f ) 


*l) ,lc-l 
1 ^ r-1 


C. + 


k.1 

r-1 ^2 ^3 


( 2 - 58 g) 

(2-5811) 


^ (1- \y 


^2 + ??T ^ 2 )'^ 


2 ,k-1 

2 ^3 


...( 2 - 581 ) 


where 


C. = 


.a. 


^1 - 2(2q+‘i-r) ’ ^2 " 2( X,- X 2 )( 2 q+ 1 -r) ' ^3 

Now 


N(i- ^2) 


= 


- X2(1- x^) 

2( Xi-X2)(^<l+''-r) 


J, Pjk’ “k = 3(p^,-p“g) (Pj2-Pj!,) * <Pj3-Pj5> (2-59) 

The different values of 3 correspond to the chain 

starting in different initial states, and, thus, different 

values of Lj^ are obtained. As mentioned earlier, in order 

2 

to obtain an 'ensemble' average the calculated 4 ® 
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weighted by the stationary state probabilities of the 
corresponding initial starting state. 

For the six steady states, probabilities are given 

bel ow 


a 2 1 - 1 ’ q 

2(2q+1 -r) ’ 2(2q+1 -r ) » 2(2q+1-r) ’ 2(2q+1-r) ’ 2(2q+1 -r) 


2(2q+1-r) 


Substituting these values, we get, 


T _ q T , ^ T 

k “ 2q+1_r ^1k + ^2k q ^3k 


(2-60) 


On substituting various values following expression for Lj 
is obtained 

>k 


k 2 q+ 1 -r 


, 2 k k 

5 (2q-1+r) + 32C, A, + 32C, 


(2-61) 


Now, 


n^1 

I X 

k=o 


n X 


114. 1 


(1-^)^ (1- 


^ 1 / -I \ ^ -t 

I (n,k) (-1)“ = A .2. hll, ±1 

k=o 


= A = 2 + 


( 2 - 62 ) 


(2-63) 


Substituting equations (2-6l - 2-63) in equation (2-56) 
and considering large values of n, the value of the expected 
square of the end-to-end distance for an oriented hydrocarbon- 
type molecule is reduced to a simple expression 

.2 


<E„> = 


* 4(2q!i?f ' ) ' ^ tU(Q-r)(1-q)j 

...(2-A) 


±2 

3^0 
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It should be possible to test the validity of this 

equation from measurements on light scattering, viscosity, 

etc., for hydrocarbons of various chain lengths. This would 

also permit an estimate of the various transition probability 

2 

parameters and from measurements of <R > on chains of 

n 

different molecular weights, the self-consistency of these 
parameters could be checked. 

The result for the special case of Tobolsky^ is 
obtained by substituting q = — , i.e., p = q. In 

this case, the expression for (equation 2-^4) reduces 

to 

= n?- 5 (2-65) 

r = 1/3 implies that there will be no energy 
difference between the trans and the two gauche configu- 
rations. Substituting r = 1/3 in equation (2-65), we get 

<R^> = 2nl^ (2-66) 

no . 

and r = o implies that the trans oonf igurations are nearly 
excluded. Substituting r = o 

<R^ = fnl2 (2-67) 

These results are well known in the theory of 


random walk. 



In this problem, the transition probability matrix 
of the Markoff chain relates the orientation of the vector 
^i ^i-1 ^i-2’ Correlations in the more distant 

links are disregarded. As a result of this limited memory 
of the links, it is not possible to prevent the occurrence 
of conformations containing loops or of conformations that 
are ruled out because of steric hindrance. This influence 
of excluded volume effects the average square length of 

alkane-type chains and this has been studied by Hijmans and 

7 

Holleman by excluding cyclohexane type loops and configu- 
rations involving pentane type steric hindrances. Such 
calculations approximate the physical situation much better. 
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CHAPTER - III 


PHASE TRANSITIONS IN POLYMERIC SYSTEMS 

3.1 INTRODUCTION 

A system is said to undergo phase transition when 
there is a sudden change in its properties with respect to '■ 
some thermodynamic parameters such as temperature, concen- 
tration, etc. A transition is said to be of the nth order 
when the nth derivative of free energy with respect to any 
of its natural thermodynamic variables has a discontinuity. 
For example, melting of crystal, helix-random coil transi- 
tions, etc. are termed as first order transition, whereas 
superconducting transitions in zero field and X transitions 
like order -disorder transformation in g brass If and 
transition in ortho hydrogen are termed as second-order 
transitions. In the present work, the problem of phase 
transition in some finite macromolecular systems is studied. 

In case of helical polypeptides, this problem has been 

1 2 

studied by Zimm and Bragg , L if son and Roig , Gibbs and 

Di Marzio"^ and others ” . The theoretical treatment is 

based on one -dimensional models. 

Actually, a polymer is a three-dimensional system 

• but it is observed that in the neighbourhood of transition 

a one -dimensional model gives a good description of its 

7 8 

properties* It has been “ shown * that a one -'dimensional 
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system does not exhibit phase transition, if it is assumed 

that the interaction has a finite range and is non-singular . 

If these two restrictions eire relaxed, then the system does 

show a transition from one to the other state. Based on 

these considerations, which is essentially a modification 

1 2 -< 

of the Ising model, Zimm and Bragg besides several others 
have given a theory for the helix-coil tramsition in simple 
polypeptides . 

The polypeptide chain consists of amide groups 

(D 

— N — C 

! 

H 

connected by intermediatry ceorbon atoms. 

Tn the a-helix, the hydrogen atom of each amide group 
forms a hydrogen bond with the oxygen atom of the third 
preceding amide group, i.e., each segment is bonded to the 
fourth subsequent segment as shown in Figure 5 * Here an 
amide group plus one adjacent carbon is referred to as a 
segment of the chain. 

Formation of a helical region is a cooperative 
phenomenon, i.e., once first segment is bonded to the 
fourth, it is relatively easy for the second segment to 
bind to the fifth and third to the sixth and so on. It 
is difficult for the first segment to bind due to a large 
loss in entropy. Once foi*tnedj however, this, in turn. 




BOTH A RIGHT* 
X ARE SHOWN 
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acts as a nucleus to which further turns can be added by- 
hydrogen bonding. This difficulty in binding is expressed 
by a parameter called the nucleation parameter a • A 
change of a few degrees in temperature or a few per cent 
in solvent composition can convert the randomly coiled 
form to the a-helix. Actually, this change in temperature 
and solvent composition brings a decrease in entropy which 
is essential for the formation of hydrogen bond. 

Zimm and Bragg^ have done the problem of transition 
between the helical and randomly coiled forms of a polypeptide 
chain. Two principal parameters introduced are (i) a sta- 
tistical one that is essentially an equilibrium constant 
for the bonding of segments to a portion of the chain 
which is already in helical form, and (ii) a nucleation 
parameter. The partition function is calculated from the 
Ising matrix as the nth power of its maximum root. Helicity 
as a function of temperature is calculated for pclybenzyl- 

glutamate and is in good agreement with experimental data. 

g 

Poland and Sheraga"^ in a recent publication have given a 

comparison of theories of helix-coil transition in poly- 

2 

peptides. Lifson and Roig have used -a sequence generating 
function whereas Pell Or ^ has used a combinatorial approach. 

In the former, a cubic secular equation is obtained, the 
maximum root of which is approximated by a quadratic 
equation. If the use of triplets eis units is replao-ed by 



the use of amino acid residues and the positive energy of 
formation of a hydrogen bond is assigned to single helical 
states, then the Feller's result is identical to the 
quadratic approximation. 

3.2 CELL MODEL FOR DIFFUSE PHASE TRANSITION 

The cell model has been applied to explain x point 
anomaly found in ammonium hexaf luro- silicate , NHj^Cl, 

KH2P0^, KHgASO^ , Ag 2 H^IO£ and . It has also been applied 
to explain thermally induced phase transition in a number 
of single strand and double strand helices^^’^V In the 
present work, the cell model has been used to study melting 
phenomenon in finite systems. This model can be briefly 
explained as follows. 

Consider a lattice having two phases, A and B, 

dispersed in each other. Assume that each region of phase 

A or of phase B is wholly enclosed by a boundary and that 

the same coordination number Z is valid for their lattices. 

Suppose that the pair partition functions for aa, ab and 

bb pairs are f^^^ and f^^ respectively. 

Consider a cell model for the tctal lattice of n 

molecules. For a molecule in the interior of one of the 

phase A regions, the contribution to the total partition 
z 

function is f ' : for a molecule on the boundary of a 

3,Si 

phase A region, the contribution to the lattice partition 

z z /h' 

fixnction is (for simple boundaries) f^^^ ' f^^ ' ; for a 



molecule in the interior of one of the phase B regions 

HZ / 2 

the contribution is f , ' and for a molecule on the 

bb 

boundary of a phase B region, the contribution is 

z/4 z/4 

bb ab 

The partition function for the lattice as a whole 
is solved by the Ising matrix method. The matrix is 


P = 


z/2 


z/4 z/4 


aa 


aa 


ab 


z/4 z/4 

bb ^ab 


z/2 


bb 


(3-1) 


The partition function Q for the entire lattice is 


Q s A* ♦ 

I 

where and X^ are the roots of matrix P. 

3 .3 MELTING OF N-PARAFFIN HYDROCARBONS 

The present modification consists of expressing 
the Ising matrix in terms of segment pair partition 
functions which are either in the amorphous region (denoted 
by r) or in the crystalline region (denoted by h) . The 
boundary between the crystalline and amorphous states, 
i.e., the first h segment in a sequence of h states is 
denoted by a k state, f ^ , fj^ and fj^ are the corresponding 
segment partition functions of the three states . Taking 


(3-z) 



into account the near est -neighbour interaction, the basic 


matrix for 

the 

problem is 






r 

k 

h 


r 

f'^ 

r 

f^' 

r 

i 4 
f2 ff 
r k 

0 


P = k 

r 

1 

k 

0 

1 

f^ 

k 

4 

h 

1 

fZ 

I h 

1 

r 

0 

4 

1 

k 


The partition function for the entire system is 
obtained and the matrix method is used to obtain fractional 
crystallinity and average sequence length in the amorphous 
and crystalline regions. The results obtained are valid 
for an arbitrary and finite value of N, where N is the 
number of segments in the chain. This enables one to 
study the melting phenomenon in normal paraffin hydro- 
carbons . 


The eigenvalues of P are given by 







= b 


(3-4) 


The matrix T which diagonal ises P is given below 
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1 

1 

T = 

N 

-f 

r 

2 

1 1 

r k 

f^ 

r 

f^ 

k 

J_ 1 

r k 


4 

1 

f^ 

r 

4- 4- 

1 1 
_fT f-2 

r h 



-^h 


Similarly 

T can be written as 





1 1 

C,f^ tl 
Irk 

^k 4 





VTv^ 

^-1 


^2 

1 1 

1 1 

^2^1 ^k 

± 


^2 

^( ^2‘^h 



0 

"3 

1 1 

P p 2 -p ’"'2' 

'^3 



N - ^h 

^2 


where 


~ \ \ 


^ - ^2 

> ^ 


and 


1 




f^ 4 

r k h 

- f"^ 

r k 



(3-5) 


(3-6) 


Let the first segment be assumed in the amorphous phase 
(gauche configuration). Its contribution to the partition 
function is given by 


U = (fj, 0. 0) 


(3-7) 


The way the matrix P is constructed, the configurations 

1 X 

ending with r, k or h have to be multiplied with f ^ , f^, 
f^ respectively before their contribution is summed to 
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give tlie partition function. For this purpose, a column 
matrix U is defined. 


U 


1 

f^ 

r 

± 




(3-8) 


Now the partition function is given by 


Z = u"^ = UTI^^ TT~^ U"*" 


( 3 - 9 ) 


From equations (3-3) to (3-8) 


UT = ( fj, f|, fj ) 


T-' 


C, f"^ 
1 r 


=2 


0 


and 


^-1 „^ N -1 
T PT 


.N-- 


0 

0 


0 


.N-1 


0 


0 

0 

0 


Therefore, the expression ( 3 - 9 ) for partition 


function becomes 



55 


UT 


[T--* PT] 


= ( f2 ) 

^ r r r ^ 


t"^ u"^ 


xN-1 

0 

0 


■^1 


‘i 

T 

0 

,N-1 

^2 

0 


^2 


1 

’r 

0 

0 

0 



0 

j 


(3- 


or 


Z = 


C. 

Z. 


r 

^1+0^ 
+ Zo 


,N 


(3-11) 


Let the number of r segments, k segments and h 

segments (exclusive of k segments) be denoted by n^ , n^^ and 
* 

nj^. These values are given as follows, 


and 


N 


2 

N 


31nZ 

9lnf 

r 
9X 


1 

Xi 


3f 


1 ^ 


3X, 


af 


NZ 


Z^ 31nC^ 
3f 


r 

Z^ 3 InC^ 
3f 



f 

k 

r^, 

3Xi Zg 

3 X 2 ' 

^k 

Z^ 3lnC^ 

BlnC^ 

rj ^ 

N 

" Z 


^^k * h. 

'^kj 

NZ 

L 

* ^2 3f, 

k J 


(3-12) 


(3-13) 


* 

"h 

^h 


”1 ^2 

3 X 3 “ 

^h 

1 31nC. 

31nC^ " 

N 

= 2 


^ ^h ^ ^2 

S^hJ 

NZ 

1^1 3 fj^ + ^2 

3f‘h_ 


(3-14) 


These expressions finally reduce to the following, if we 


10 ) 


substitute 
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and 


^h 

f - Y 

(3-15) 

r 


^k 

f - ° 

( 3 - 16 ) 

r 


A = Zyz 

(3-17) 


( y P 


(3-18) 


= ^ + (2A-1) + [(2A-1)p -I+Y] 

2 P'^N 


(3-19) 


N 

■H- 

N 


- (y (1+Y+ p ) A a (1+Y) ( 1-y) / ^ o^^ 

2p!y- o) ^ 1y- =>)P ■ 2h'p^ C(2A-1)P-1+Y ] (3-20) 


I [ l-Y-t-P .1 

P 1+Y-P 


AY(l-Y-f2 o ) 
p(Y- a) 


NP 


[(2A-1 )P-1+y3 


- . .(3-21) 


It can be immediately checked that 


n 


N 


k 


N 


* 

N 


= 1 


( 3 - 22 ) 


Inasmuch as we use segment peirtition functions, 
these quantities can be interpreted in the following manner, 


II 

(3-23) 

AHf 


ffa = «h < RT ) 

(3-24) 


(3-25) 
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where is the statistical weight of an r segment, gj^ 

is the statistical weight of an h segment and A is the 

heat of fusion from the crystalline to amorphous state. 

The partition function fj^ of the k segments which 
represent the boundary between amorphous and crystalline 
regions is taken to be temperature independent and very 
small so as to represent the difficulty encountered in 
entering such a boundary region. In equations (3-23 - 
3 - 25 ) > gj^, and e are assumed to be constants, and Ah^ 
is the heat of fusion per mole of segments. For large N 
at the fusion point, fj^ is equal to f ^ . The temperature 
of fusion and the entropy of fusion are, therefore, given 
by 


wher© 


Hence 


where 


and 


T^ = AH^/Rln(gygj^) = AH^/ A (3-26) 

A = Rln (gj,/gj^) 

Y = ^ exp( AH^/RT) (3-2?) 

r 

= exp [( AH^ - AS^T)At] 

AH„ „ 

= 8XP [ 55 ^ (T^ - T) ] = e* (3-28) 


X = 


AH^ 

TT^ 

f. 


(T„ - T) 




The number of amorphous sequences r and that of 


helical sequences h are given by 
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4 ^ 


(3-29) 


(3-30) 


Let Q be the fractional crystallinity, then 


Q = 1- 


(3-31) 


= (1-Q) 


(3-32) 


Now the density d of a macr omolecule with large N 


can be written as 


QO + (1 - Q) A 


(3-33) 


where C and A are partial densities of crystalline and 
amorphous parts . 

Value of C can be obtained by taking the crystalli- 

12 

nity and density values at N = 30, as the experimental 
density curve is available for this value of chain length 


Now at 


T = 273°k 

Q = 0.865 

d = 0 , 952 . 


This gives 


c = 0.979. 


At sufficiently high temperatures. 


Q = 0 


therefore, A = d 
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where is the density of the amorphous part and its 

1 3 

value as given by Flory, Orwoll and Vrij is 0.779 • 
Substituting the values of C and A in equation 
(3-33), we get 


d = 0.2Q + 0.779 (3-3^) 

Normal paraffins for odd N ^9 and for even 
N > 44 exist in the orthorhombic crystalline form. Transi- 
tion from the orthorhombic to the hexagonal form precedes 
melting for odd N in the range 1 1 IN £.43 and for even N, 

34 <N £.42 . Flory and Vrij^^ have reported the melting 
points £ind enthalpies of fusion for N = 15, 19, 25, 29, 30 
and 43 . Their data is given in Table 1. 


- TABLE 1 


Chain length 

N 

Enthalpies of 
fusion Ah 

Melting point 
ok 

15 

700 

280.5 

19 

750 

302.5 

25 

800 j 

324.5 

29 

805 

330.0 

30 

795 

336.0 

43 

800 

358.0 


Use has been made of the above data in obtaining the 
melting curves Q versus T. The calculated data is given 


in Table 2 . 




TABLE 2 

CRYSTALLINITY VS. TEMPERATURE FOR 
VARIOUS CHAIN LENGTHS AND cr VALUES 


Temp-e- 
ratur e 

Cry 3 1 al 1 ini ty 

t - - 

°IC 

N = 15 

N = 19 


a =10"^ 

CT = 2 x 10 “^ 

a= 10 "^ 

,^-2 
cr = 1 0 

0 = 6 x 1 0 "^ 

a= 10 -'' 

250 

0.747 

0.717 

0.348 

0.814 

0.812 

0.769 

260 

0.719 

0.669 

0.241 

0.795 

0.783 

0.694 

270 

0.687 

0.611 

0.163 

0.773 

0.742 

0.587 

280 

0.652 

0.546 

0.110 

0.748 

0.686 

0,457 

290 

0.613 

0.477 

0.075 

0.720 

0.612 

0.330 

300 

0.572 

0 . 4 o 9 

0.052 

0.688 

0.524 

0.226 

310 

0.529 

0.346 

0.037 

0.652 

0.430 

0.151 

320 

0.487 

0.290 

0.027 

0.612 

0.341 

0.100 

330 

0.445 

0.242 

0.020 

0.570 

0.263 

0.068 

3^0 

0.405 

0.202 

0.016 

0.527 

0.201 

o.o 47 

350 

0,369 

0.170 

0.012 

0.483 

0.153 

0.033 

360 

0.335 

0.143 

0.009 

0.441 

0.118 

0.024 

370 

0.304 

0.122 

0.008 

0 . 4 oo 

0^091 

0.018 

380 

0.277 

0.105 

0.006 

0.362 

0.072 

0.014 

390 

0.253 

0.091 

- 

: 0.328 

0.057 

0.010 

4 00 

0.231 

0.079 

- 

0.297 

0 .o 47 

0.008 

4 io 

0.212 

0.070 

- 

0.269 

0.038 

0.007 

420 

0.196 

0.062 

- 

0.245 

0.032 


430 

0.181 

0.056 


0.223 

0.027 

- 

44 o 

0.168 

0.050 


0.205. 

0.024 ‘ 

- 

450 

0.157 

o.o 46 


0.188 

0.020 

- 

46 o 

0.147 

0.042 

- 

0.174 

0.018 

- 

470 

0 

00 

0.039 

- 

0.161 

L 

0.016 

- 


( c ontd . ) 














Table 2 (contd.) 


Teti 5 )e. 

r atur < 

3 


Crystallinity 



°K. 


N = 30 


1 

I 

N = 43 



a =10“^ 

-5 -5 

a = 9 x 10 '^ 0=10 

j - - 

a = 1 o~^ 

CT = 5 x 10 “^ 

cJ =10“^ 

250 

0.892 

0.893 

0.893 

0.925 

0.926 

0.926 

260 

0.881 

0 ♦ 882 

0.881 

0.917 

0.918 

0.918 

270 

0.869 

0.869 

0.864 

0.908 

0.910 

0.910 

280 

0.854 

0.854 

0.835 

0.898 

0.900 

0.900 

290 

0.838 

0.832 

0.774 

0.887 

0.889 

0.889 

300 

0.817 

0.798 

0.654 

0.873 

0.876 

0.874 

310 

0.791 

0.742 

0.471 

0.857 

0.859 

0.853 

320 

0.758 

0.648 

0.283 

0.837 

0.835 

O.813 

330 

0.712 

0.518 

0.150 

0.812 

0.797 

0.722 . 

3^0 

0.652 

0.389 

0.076 

0.781 

0.726 

0.547 , 

350 

0.575 

0.240 

0.039 

0.739 

0.597 

0.325 ' 

360 

0.486 

0 . 149 

0.021 

0,682 

0.418 

0.158 

370 

0.394 

0.091 

0.012 

0.606 

0.248 

0.071 

380 

0.308 

0.057 

0.007 

0.512 

0.133 

0.033 

390 

0.236 

0.034 

- 

o, 4 o 9 

0.071 

0.016 

4 oo 

0.179 

0.025 

- 

0.312 

0.038 

0.008 

hio 

0.137 

0.017 

- 

0.231 

0.022 

- 

420 

0.106 

0.013 


0 

• 

0 

0.014 

f 

430 

0.084 

0.009 

- 

0.126 

0.009 

- 

44 o 

0.067 

0.007 

- 

0.096 

- 

- 

450 

0.055 

- 

- 

0.075 

- 

- 

46 o 

0.046 


- 

0.059 

- 

- 

470 

0.039 

- 

- 

o.o 49 


wm 
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Tlie nucleation parameter is adjusted to give 
correct melting point. Crystallinity versus temperature 
curves are shown in Figures 6a-6f . The dependence of 
transition cmrves on cf values is also shown and, as expected, 

it becomes sharper for small values. 

1 2 

Seyer , Peterson and Keays have reported the 
variation of density with temperature. These results can 
be compared with the values of density obtained from 
equation (3-3^) • The experimental results and theoretical 
calculations for N = 30 are shown in Figure 7- The corres- 
ponding data is given below in Table 3» 


TABLE 3 



Because of the equivalence of the crystal structures 
for polyethylene and N-paraffin hydrocarbons and the virtual 
identity of their unit cell dimensions, it is expected that 
the limiting melting point for N-paraffin hydrocarbon aeries 
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must be coincident with the melting point for linear 
polyethylene of infinite chain length. This is supported 
by the curve shown in Figure 8 . The value of the nucleation 
parameter also asymptotically approaches the value 5x10“ 
(Figure 9) i which yields correct results for polyethylene. 
Table k gives the values of nucleation parameter for 
different chain lengths. 


TABLE 4 


Chain 

length 

N 

Nucleation 

parameter 

15 

_3 

2 X 10 

19 

1 

o 

25 

-4 

2 X 10 


_4 

29 

1 X 10 

30 

9 X 10”^ 

j ^3 

5 X 10“^ 


3,4 EFFECT OF PRESSURE ON TRANSITION IN POLYETHYLENE 

The mechanical behaviour and morphological stiiicture 
of a polymeric material are markedly changed when subjected 
to stress. For example, the elongation of a network polymer 
in the amorphous state increases its free energy and if the 
polymer is capable of crystallinity, its melting^ point will 




MELTING POINT (in°k) 




(>9 


also be increased. Vulcanised rubber which, does not 
crystallize at ordinary temperature becomes highly crystalline 

when sufficiently elongated, even at temperature above 100°C. 

1 5 

Matsuoka has argued that the overall specific volume below 
the melting point is not a unique function of pressure and 
temperature and hence may not be regarded as a thermodynamic 
quantity. However, a temperatur e -pres stir e relationship is 
very useful in estimating the true melting point as a function 
of pressure. The problem regarding the effect of pressure on 
melting point has also been studied by Krigbaum . The 
melting point increases as the applied pressure is increased 
and can be represented by equating Gibb's free energy of 
fusion to zero. For segments within a sufficiently large 
crystalline domain, 


AG = AU + P A - T A S. 


u 


(3-35) 


where aG is the free energy, A U is the internal energy 
of fusion, P is the pressure, AV^ is the volume of 
fusion, T is the absolute temperature and A is the 

entropy of fusion. The change in melting temperature due 
to the applied pressure is given by Clapeyron's equation 


dT^/dP = T^ AV^ / ah, 


u 


(3-36) 


where AV and AH are the change in volume and change 

*ul 



70 


in enthalpy at the melting temperature T . 

m 

In the present work the formalism of the X&xng model 

for a linear chain of spins in a magnetic field has been 

study the effect of pressure on the transition 

curve for polyethylene* When a polymer segment is subjected 

to stress, the contribution of a segment to partition function 

in crystalline or amorphous region is also altered* In the 

absence of any pressure, the contribution of the segments 

in the disordered and ordered state are f and f, res- 

r n 

pGctively and for the first ordered state in a sequence of 
ordered states, this contribution is f^^. When pressure is 
applied, the contribution of a segment in the crystalline 
or helical state relative to an amorphous state may be 
written as ffj^ where f = e“ being the energy 

difference between the two states due to. the pressure alone. 
The corresponding Ising matrix now becomes 


r 

yr = k 
h 


f f, 0 

r k 


0 

0 


ff,. 


ff 


h 


(3-37) 


The partition function for the entire system is obtained 
by the same procedure as done in the previous problem. 

The eigenvalues are given by the following relations 
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(3-38) 


X. X = ff f, - f f 
12 r h r k 


(3-39) 


The matrices T and T” are 


(3-4o) 




1 h 




(3-41) 






C.ff. f, 

1 h k 


^1 


C^ff, f, 
2 h k 


^( 


(3-42) 


T" T = 1 gives 

C - C 


^2-^^h 


, C, 


X,- X^ ’ -3 


h. k 




...(3-43) 


Initial and final states of the system are given by 


U = (f^, 0, 0) and U'' 


(3-44) 
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because we consider that the initial segment can be only in 
the random state and the final segment can be in either of 
the three states . 

Partition function of the system is given by 


Z 


Utt-I u’’’ 


XN 

1 


XN 

2 


(3-^5) 


Let the fraction of r segments be denoted by 
then 


< n > 
r 


1 9lnZ 
f 


^ a Z 

N 9lnf “ NZ ^ 
r r 


NZ 


VO tlN-1 ^ ^1 ,N-1 


9X 

TT 


2 ,N 'S ""^2 

+ IT" 2 a f _ 


c x^-"' 


L-^o \N-1 2 

9f 2 9f 

r r 


JSZ 


iN 


9 X^ 9 X 2 

Substituting the values of — and -g £- , 


=’<=1 «2 

‘l 3f * 2 8f 

r r 

...(3-46) 


<11 > 
r 


^ ^ 1 " 


.ff _ ^Jll) 

■1 “h ^ xj 




1 + 


(^) 


N 


NZ 


9C 

N .N 

9f ^2 


a c. 


TT 


r JL 


( 3 - 47 ) 


For N large, the above expression reduces to 


<n^> 


^ - ^^h 

( N - 


(3-4fi) 
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Degree of order or crystallinity is given by 


Q 


1 - < n > = 

r 


- Xg + ff^ 

Xi - X2 


(3-49) 


Substituting 


— = o, 




^ _ V., ^ = Y2, = yi 

r ^r ^r 


Sind 


"k ■ 


= <3 


equation (3-49) becomes 

Q 


where YI 


and 


Also 


fS = 


fS - Y2 
YI - Y2 


= j j^1 + fS + *^(1-15)^ + 4a S 


Y2 = 


2 

2 


S =? exp 


1 + fS - (l-fS)^ + 4cr s 


■ah(t^-t)' 


RTT^ 


, AEv 
exp(— ) 


( ah + AE)/RT - AH/RTf 
© © 


(3-50) 


AH'/iiT - AH/RTf 

zz 0 © 


( 3 - 51 ) 


where AH' = AH + AE 


Degree of order as a function of temperature at 
different values of pressure is calculated from the expre- 
ssion (3-50). The heat of trsinsition per mole for a residue 
from the ordered to the disordered state (trans to gauche) 



at one atmospheric pressure is 943.0 Cal. /mole. The 
transition temperature is 4lO°K. The entropy change 
associated with the transition is 2,30 Cal./mole/oK ar i d 
it is assumed to be the same at all pressures. This is 
justified from the experimental data on entropy of fusion 
as a function of pressure reported by Matsuoka^^. The 
deviations from the average value of 2,30 are of the 
order of 0.10, The heat of fusion, AH' in generating 
transition curves at different pressure are given in 
Table 5. These values are taken from experimental results 
reported in reference 15 . 


TABLE 5 


pressure 
(in atmospheres) 

Heat of fusion 

AH' 

(Cal ./mole) 

Melting 

point 

(°K) 

1 

943 

4lo 

1 130 

952 

4l4 

1 1 000 

1005 

437 

I 2000 

1053 

458 

1 3000 

i 

1101 j 

1 

479 


The calculated values of degree of order vs, 
temperature are given in Table 6 and represented in 
Figure 10. The general characteristics of melting under 
the atmospheric pressure are also observed at elevated 




TABLJ5 6 

DEGREE OF ORDER VS. TEMPERATURE AT ViiRIOUS PRESSURES 



0.999 

0 . 99 '^ 

0.977 

0.500 

0.024 

o.oo6 

0,003 

0.002 


0.999 

0.996 

0.992 

0.963 

0.218 

0.016 

0.012 

0.002 


0.995 

0.987 

0.884 

0.059 

0.010 

o.oo4 

0,002 


0.999 

0.999 

0.998 

0.997 

0.996 

0.989 

0.931 

0.1 04 

0.013 

0.005 

0.002 


0.999 

0,998 

0.998 

0.996 

0.990 

0.95^ 

0.196 

0.018 

0.006 

0.002 








temperatures. The shift in melting point with pressure 

is given in Table 5. It is observed that the shift in 

tlie melting point with increasing pressure is exactly 

1 5 

the same as reported experimentally * The general features 
of the transition curve resemble a step function. It is 
difficult to make a quantitative comparison with the 
specific volume versus temperature measurements of Matsuota 
at different pressures because the dependence of specific 
volumes for the ordered and disordered state on pressure 
cannot be easily taken care of. This is partly because 
the amorphous state itself cannot be uniquely defined and 
it is evident from the fact that Matsuoka’ s measurements 
of the specific volume of the melt varied considerably 
from sample to sample. Further, it is nowhere specified 
by the author whether his measurements were carried on a 
well annealed sample or not, A detailed comparison with 
his transition curves, therefore, becomes difficult because 
the history of the sample is not known. An initial distri- 
bution of the crystallite size would tend to broaden the 
trsuasition. Branching has a very large effect on the 
sharpness of the transition. Due to this the transition 
becomes very diffuse. 
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3.5 HELIX-COIL TRANSITION IN A TRICOPOLYMERIC SYSTEM 
The cell ciodel for diffuse phase transitions 
described previously can be modified to explain helix-coil 
transitions in two and thr ee-copolymer ic systems also such 
as (glutamic acid-lysin)^ and (glutamic acid-lysin-alanin-)^. 
The modification consists in introducing two or three random 
and helical states respectively. Considering it specifica- 
lly for a thr ee-copolymer ic system, let the three helical 
and random coil states be designated by h^.h^.h^ and r^, r^ , 
r^ and the corresponding segment pair partition functions 
by f, , f, , f, and f , f , f respectively. The k 

hi h^ h^ ri Tg r^ 

state which represents essentially the boundary tension 

between helical and random coil states remains the same. 

Any difference arising due to the different helical states 

is taken care of by the nucleation parameters — — , ^ 

fr k k 

and — — . The Ising matrix for the partition function 

can be written as 
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"’2 

"3 


^2 

iij 

k 

"’l 

1 1 
f2' f2 

i 1 

f2 £2 

T r 

12 

1 1 
£2 £2 

T r 

13 

0 

0 

0 

ft 

r ^ k 

^2 

4- J- 

fT f2 
^2 

1 1 
f2 f2 

^2 ^2 

1 1 

f? f2 

2? r 

2 3 

0 

0 

0 

1 1 
f2 f2 

^2 ^ 

^3 

1 1 
f2 fS 

^3 

f’^ f"^ 

^3 ^2 

1 1 
f.2 f2 

r3 

0 

0 

0 

1 1 

f2 £2 

k 

^1 

fi fi 

N ’'i 

1 1 
f 2 f2 

^1 ^2 

1 1 

h^ 

1 1 

hi h^ 

1 1 
f2 f^ 
hi h^ 

1 1 

il. ^ xl.^ 

0 

^2 

1 1 
f? £2 

^2 

1 1 

f; 

hj 

1 1 

^2 "^3 

h2 1^1 

fi fi 
h2 h^ 

1 1 

^2 ^3 

0 

^3 

f i f i 
h3 r^ 

1 1 
fj f2 

^3 ^2 

1 1 
f 2 fT 

hj r3 

1 1 
f^ f2 
h3 h, 

1 1 

^3 ^2 

1 1 
h^ h^ 

0 

k 

fl f ^ 
k r^ 

ft f 

k rj 

ft f* 

k 

f^ f^ 
k h^ 

1 1 
f^ f'^ 

^ ^2 

1 1 

f J 

k h^ 

0 


...(3-52) 


The matrix M has five zero roots and of the two remaining j 
the greater eigenvalue is given by 


max 


= r *‘'r *^h +^h > + ^ +‘'li 

L ^2 3 1 2 3 1 2 3 1 2 3 


Define 


-4 

-..(3-53) 


f + f + f = 
^2 ^3 


(3-54) 


and 




(3-55) 



sc 


Ther ef ore , 


'max 


2 

2 


+ ‘‘Vr ] (3-5«) 


The expression 0 for the fraction of segments in helical 
states is 


X 


0 = 1 


1 - + 2 + A 


A ( 1 + e + A) 


where 


(3-57) 


A = 


(e^ - 1)^ + 


X = ( Ah^ _ T AS^ )/RT, 


AH^ and A being respectively the standard heat 
of entropy of transition per mole of segments from helical 
to random coil states, 

R is the gas constant, 

T is the absolute ten^erature ; 

AH^ 

at the transition temperature = “Xs” (3-58) 

17 

In a recent publication Morita, Simon and Blout 
have studied transition in (glutamyl-lysyl-alanyl)^ trico- 
polymer by measuring optical rotatory dispersion. It was 
observed that at PH3» the copolypeptides containing 20 mole ^ 
or more alanine exhibit a helix content equal to the sum 
of their alanyl and glutamyl residue contents whereas at 
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PH8 the helix content equals the alanyl content when the 

latter is 4o mole ^ or higher. At lower alanyl contents the 

electrostatic interaction between charged glutamyl and lysyl 

residues makes some contribution to helicity. There is helix 

stabilisation as a result of electrostatic interaction^ In 

addition, hydrophobis interaction also contributes to helical 

stability. At higher PH value of 12 the amount of helix 

content observed is proportional to the mole ratio of alanine 

residues present in the polymer . In fact , at all stages the 

total helical content is the sum of contributions from alanyl, 

glutamyl and lysyl residues. The alanyl residues are helix 

promoters at all PH values whereas the contribution of the 

other residues depends on their state of ionisation and 

17 

separation. Morita, Simon and Blout have studied a tra 

copolymer with 1:1:3 mole ratio of glut amyl -lysyl— alanyl 

/ , -. 20 , 20 60. 

residues in 8M urea and 5*5 LiBr solution (poly glu ly-^ ala j . 
When the solution is heated over a temperature range from 
5° to 80°C , it undergoes a diffuse transition. At PH = 3> 
the melting temperature is about 52°C and 39®C at PH = 8. 

This gives an overall enthalpy change Ah of approximately 
-6.4 K Cal. /mole in both the cases. This enthalpy change 
can be calculated from the slope of 0 versus l/T. The 
near equality of Ah for the transition curves at PH3 
and PH8 shows that they have identical slopes. The 

AH^, associated with the 


standard enthalpy change. 
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transfer* of a single residue from random to ordered state 
is calculated from the knowledge of nude at ion parameter 
(f'kAr)- This, in turn, depends on the system and the 
solvent . In our calculations 

1 

AH^ = AH ( ^ )2 (3.59) 

r 

3 

with jr- = 3 X 10"'^ ' 

r 

and Ah = -6.4 K Cal. /mole 

is calculated from equation (3-58) by using 

= 520 c at PH = 3 

and T = 39°C at PH = 8 . 

m 

Using these thermodynamic parameters and expression 

(3-58) helicity is evaluated as a function of temperature. 

The results obtained are shown in Figure 11. They are 

in good agreement with the experimental results of Morita, 

1 7 

Simon and Blout . The circles and triangles represent the 
experimental data of these workers for PH3 and PH8 values 
respectively. The experimental results and theoretical 
calculations are given in Table 7* The calculated curve 
shows slight deviations from the experimental points below 
295°K. This could be because of the use of approximate 
and extrapolated data in using the same value of AH for 


both PH values. 
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CHAPTER - rV 


CONCLUSION 

In the preceding chapters two main problems on 
configurational statistics and phase transition have been 
discussed in two different types of polymeric systems. 
Both theoretical approachps are subject to several appro- 
ximations and hence the results obtained have a somewhat 

limited validity. Taking the first problem regarding the 

2 

calculation of expected square of end-to-end 

distance, the limitations arise from 

(i) the model considered 

(ii) exclusion of excluded volume effects 

(iii) ignoring of chain foldings. 

Two models, based on random walk on cubic and 
diamond lattices have been considered. Of the two models, 
diamond lattice model is a better approximation because 
the bond angle (109°28') is more close to the correspond- 
ing angle in polyethylene and hydrocarbons than the bond 
angle (90°) in a cubic lattice. This makes the descrip- 
tion of crystalline and amorphous regions by random 
walks with and without a preferred scheme all the more 
approximate on a cubic lattice . The situation is further 
complicated because the amorphous state itself cannot be 
uniquely defined. The use of other lattices which are 



likely to simulate the polymer structure better would 
complicate the problem mathematically. 

The second approximation is an important one. The 
interaction between various segments has been completely 
omitted. In particular, the fact that the same volume 
element cannot be occupied simultaneously by two units 
of a polymer chain is not taken into account. In other 
words, we have not accounted for the forces of repulsion 
which act between the different segments of the chain 

that have accidently approached each ether to a very 

1 . . 2 
close distance. Montroll , and Hijmans and Hclleman 

have, however, studied the problem in case of a square 

and a diamond lattice respectively* In the former, 

first order overlaps resulting from foiir bonds and in 

latter enclosures of the pentane type rings have been 

excluded. Better results have been obtained by a 

numerical solution using Monte-Carlo method and computer 

34^ 5 

simulation of polymer chain conf igur ation * , hdwards 

and others have also obtained analytic expressions for 
the excluded volume problem. The concept of the forma- 
tion and existence of Lamellar crystallites in bulk 
crystallized homopolymers, appears to be well founded 
and supported by electron diffraction and electron 
microscope observations^. On this basis, it is asserted 
that within an individual crystallite the polymer chains 
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are folded regularly. The structure of a semi-crystalline 

regarded the same as that of a single crystal 

with the amorphous part representing defects in the 

crystal. The models presented for the polymer structure 

in this work do not take care of these modern ideas. 

However , the theory can be modified to take into account 

these ideas of chain folding. In such a case when the 

crystalline regions are regular helices, there can be h 

bond vectors involved in the helix before axi advance of 

one unit in the crystal direction. The chain folding 

for a square lattice model has been considered by Gupta 

7 

and Gupta , Lastly, to make an effective comparison 
with experimental data, a knowledge of the distribution 
of end-to-end distance and its average value is needed. 
Besides there are complicating factors like branching. 

Coming to the problem of phase transition, the 

20 20 

results obtained for the tricopolymer of ^lut -lys 
ala^^ are in much better shape as compared with those 
for n-paraffin hydrocarbons and for polyethylene. In 
both cases it is assumed that there exist only two 
states- helical (crystalline) and random coil (amorphous) 
whereas in between these two there can be a variety of 
states. In case of polypeptides one of these states 
corresponds to the presence of a hydrogen bond (in a - 
helix) between two given segments of the polymer and the 



other' to the absence of such a bond between the same two 
segments. In order to choose an interaction between these 
two states which corresponds to the physical picture of 
the process of helix formation, it is clearly necessary 
to introduce very large repulsions. This is because the 
existence of one hydrogen bonded state makes it very easy 
for the next segment to bind itself. In other words, the 
presence of a bonded state strongly » repels* the existence 
of a non-bonded state. This is the physical significance 
of the nucleation parameter. The possibility of hydrogen- 
bond formation is due to the relative stiffness of the 
chain backbone with respect to bending modes. This gives 
a natural curvature to the polymer so that once a segment 
is bonded, the probability that the following segment will 
also bo in a position to form a bond is greatly enhanced. 
Thus, in case of a polypeptide chain, it is easy to under- 
stand the physical significance of cooperativity factors 
as well as the boundatry tension term or the nucleation 
parameter. The structural characteristics in n-paraffin 
hydrocarbons and polyethylene do not make the physical 
meaning of cooperativity factor and boundary parameter 
easy to understand. Cooperativity can be regarded as a 
result of both intermolecular and intramolecular inter- 
actions and the latter should be related to the height 
of potential barrier in going from one conformation of 
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bond to another. In case of a finite system the end effects 
also become important . In such polymeric systems the dis- 
tinction between ordered and disordered states gets further 
diffused. There undoubtedly exists a substantial proportion 
of domains that cannot distinctly be classified as either 
crystalline or amorphous. The nature of the noncrystalline 
regions is further complicated by the presence of voids as 

a direct result of the volume decrease during the late stage 

8 

of crystallization . Unlike the free volume related to 
density fluctuations in liquids and amorphous polymers, 
these voids, or the regions of zero density, cannot be 
treated as a thermodynamic quantity. 

In the case of branched polyethylene, it is fxirther 
difficult to decide how many carbons in the neighbourhood 
of either branch points or methyl groups can actually be 
counted as not participating in the crystalline lattice 
structure. In studying the effect of pressure on melting 
ciirve of polyethylene, there is yet another factor. The 
overall specific volume below the melting point is not a 
unique function of pressure and temperature and cannot,, 
therefore, be treated as a true thermodynamic quantity. 

In the present calculations, it is the variation of the 
degree of crystallinity as a function of temperattire at 
various pressures which is used. The results should be 
compared with the measurements on a linear polyethylene 



sample which has been well annealed. The history of the 
sample thus becomes an unimportant fact. In branched 
polymers, the width of transition is broadened. 

These are thus the unavoidable limitations in the 
theoretical approach x^resented in the preceding chapters 
in studying the configurational statistics and transition 
in polymeric chains. What is needed is a theory which 
can take care of not two but many more 'particle states' 
lying between highly crystalline and completely amorphous 


states . 
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